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Resonant activation of a synthetic antiferromagnet (SAF) is known to result in a dynamic running 
state, where the SAF's symmetric spin- flop pair continuously rotates between the two antiparallel 
ground states of the system, with the two magnetic moments in-phase in the so-called acoustical spin- 
resonance mode. The symmetry of an ideal SAF does not allow, however, to deterministically select 
a particular ground state using a resonant excitation. In this work, we study asymmetric SAF's, 
or synthetic ferrimagnets (SFi), in which the two magnetic particles are different in thickness or 
are biased asymmetrically with an external field. We show how the magnetic phase space of the 
system can be reversibly tuned, post-fabrication, between the antiferro- and ferri-magnetic behavior 
by exploiting these two asymmetry parameters and applying a uniform external field. We observe a 
splitting of the optical spin-resonance for the two ground states of the SFi system, with a frequency 
spacing that can be controlled by a quasistatic uniform external field. We demonstrate how the 
tunable magnetic asymmetry in SFi allows to deterministically select a particular ground state 
using the splitting of the optical spin-resonance. These results offer a new way of controlling the 
magnetic state of a spin-fiop bilayer, currently used in such large scale applications as magnetic 
memory. 

PACS numbers: 



I. INTRODUCTION 

Magnetic memory technology employs the effects giant 
magnetoresistance^ ^ and tunneling magnetoresistance^^ 
in spin- valves type sensors^. In these, the switching of the 
magnetic state is achieved by using either external mag- 
netic field^ ^ or spin-transfer-torqu4^^^^ writing. In this 
work we demonstrate a reliable way to switch a magnetic 
nanodevice employing microwave excitation in resonance 
with the spin eigen-modes in the device, which act to 
amplify the action of the microwave field. 

Synthetic antiferromagnets consist of two identical 
dipole-coupled magnetic layers separated by a thin 
spacer. The dipole coupling results in an antiparallel 
(AP) alignment of the magnetic moments in the ground 
state, which makes them similar to classical atomic anti- 
ferromagnets but with a tunable coupling strength be- 
tween the magnetic moments set in fabrication by a 
suitable choice of the material parameters and geome- 
try. When the lateral shape is elliptical, the system has 
two AP ground states, in which the moments are aligned 
along the easy-axis (EA) of the ellipse. Due to the sym- 
metry, the energy of the two AP ground states is degen- 
erate. 

The quasistaticfm^^ and dynamic^l^^ behavior of 
SAFs has been studied extensively. Due to the nearly 
perfect flux closure in the AP state, SAFs are charac- 
terized by high stability against thermal agitation^^. It 
was recently shown that the spin dynamics contains col- 
lective acoustical and optical resonance modes^^. It has 
also been shown that excitation at the optical resonance 
frequency can result in switching between the two AP 
ground states^, typically resulting in multiple switching 



or a spin-running state, since the resonant switching in 
an ideal SAF is symmetric for the two AP ground states 
so a speciflc state cannot be selected deterministically. 



In order to enable use of resonant excitation for select- 
ing a speciflc ground state, which is highly desirable for 
applications, the spin-resonance spectrum can be tuned 
such that, for example, resonant switching occurs at dif- 
ferent frequencies for the two ground states. For this 
the spin-flop bilayer must have suitable magnetic asym- 
metry, such as asymmetry in the the magnitude of the 
two magnetic moments or asymmetric fleld biasing of 
the layers externally so each layer experiences a differ- 
ent Zeeman field. The former can be achieved by making 
the two layers of the same material but slightly differ- 
ent thickness or of the same thickness but of materials 
with different saturation magnetization. The latter can 
be built-in into the nominally flux-closed reference layer 
of the read out junction^^. These two magnetic asymme- 
try contributions can also be combined, which we do in 
this study for achieving greater flexibility in controlling 
the energetics and spin dynamics in the SFi system. 

In this paper we develop a method of deducing the in- 
dividual magnetic asymmetry contributions in SFi, post- 
fabrication, and show how they can be used to reversibly 
tune the magnetic behavior of the system between SAF- 
like and SFi-like. We further show how an external uni- 
form magnetic field can be used to control the symmetry 
of the magnetic double-well potential as well as the split- 
ting of the optical spin resonance in SFi with a thickness 
imbalance and asymmetric local field-bias. 
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FIG. 1: (Color online) The top- left panel illustrates the layout of the nanopillar, which integrates an SFi bilayer and a nearly 
flux-closed read-out layer, separated by an Al-0 tunnel junction, subject to a quasistatic magnetic field Hx directed along the 
easy axis; the top-view schematic shows the magnetic layout of the bilayer, with the two dipole-coupled moments at angles Lpi 
and Lp2- The optical spin resonance for the two ground states is illustrated in the bottom- left panel for a symmetric bilayer 
(SAP; single resonance peak) and an asymmetric bilayer (SFi; resonance doublet). The magnetic asymmetry is achieved by a 
thickness imbalance (ti 7^ t2) and/or asymmetric fringing field from the reference layer acting on the individual layers in the 
spin- flop bilayer (if]^ 7^ H^). The top-right panels illustrate the symmetric and asymmetric magnetic potentials, generally 
corresponding to SAP and SFi. Interestingly, as shown in detail analytically in the text below, suitable magnetic asymmetry 
combined with a uniform external biasing field (Hx) can be used to reversibly tune the system between SAP- like (symmetric 
potential minima and a single resonance) and SFi-like (asymmetric potential and a resonance doublet) behavior. Hx — Hgf is 
the spin-flop field of the bilayer, at which one of the minima vanishes completely and the two macrospins switch into a scissor 
state. The degeneracy in the magnetic potential and the sign and magnitude of the resonance splitting depend on the sign and 
magnitude of the external uniform biasing field (see text for details), which provides an easy and powerful tool for mapping 
out and controlling the static and dynamic properties of the system. 



II. SFI EIGEN-MODES 

In this section we develop specific predictions for spin 
dynamics in SFi, which are then used to guide the exper- 
iment. In particular, we find analytically that the optical 
spin resonance in the system splits into a doublet, which 
allows controlled selection of the individual ground states 
using resonant microwave pulses modulated in frequency. 
We start by illustrating in Fig. [l]the energetics of SFi, 
including the key symmetry criteria for SAF-like versus 
SFi-like magnetic behavior. 

In the macrospin approximation, appropriate for our 
sub-micrometer junctions, with the z-axis out-of-plane 
and the in-plane easy-axis (EA) due to shape anisotropy 
along X, the total magnetic energy (E) of the bilayer can 
be written as^ 

E ^ 

+ 7jxmixm2x + ljymiym2y + 7^2^12^22 (1) 
-2{hj: + h'li.)mix - 2hy cos {ujt)miy}, 

where Ms is the saturation magnetization of the material 
(here taken to be same in both layers); Vi = Tvabti/A the 



volume of layer i with a, b and ti the length, width, and 
thickness of that layer, respectively; Ni^ the demagnetiz- 
ing factor of layer i in the a-direction, with Nia = 1 
and Ni^j^^yj = n^x^yyti/b, and the reduced demag- 

netizing factors^; the cartesian unit magnetization 
vector of the i-th layer; hu = Hu/^'nMs the normal- 
ized intrinsic uniaxial anisotropy field along the easy- axis; 
hy = Hy/AiiMs the normalized amplitude of the applied 
microwave field in the y-direction; uj the frequency of the 
microwave field; hx = Hx/^ttMs the normalized external 
quasistatic easy-axis field; = H-^/AttMs the normal- 
ized asymmetric biasing field (fringing field) acting on 
layer i in the x-direction; = r^Nia the inter-layer ex- 
change field from layer i acting on the other layer in the 
a-direction^. In our case is solely of dipolar nature. 
For simplifying the derivation to follow it is convenient 
to rewrite all terms proportional to ti through 

Q^=[l- i-iye] Q, e = {ti- t2)/{ti + ^2), 

such that Qi Nix,Niy,^ix,^iy,Vi. 

The natural coordinate system for SFi is polar, in 
which 

= {cosipi sin^^, sin(/9^ sin^^, cos^^). 
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For thin magnetic particles elongated in-plane in x, the 
following holds: 

Ni:, < Niy < Ni, « 1; ei= n/2 + < 1, 

I 

W -Wq 



Here I^o = ^[(1 + e^){Ny + TV,) + /,JM|; x 
are the normal coordinates and are given by 
^ = ((/:)! + ip2)/2, X = (^1 - ^2)/2; = ^1 + 6, 
Iz = — ^2' Further, the following notations were used: 
hf = {h^ + h^)/2, hd = {hf - h^)/2. 



and cos(^^) ^ — Quantities r^, ^i, ^2 are small param- 
eters in the problem, therefore ^zii^j can be neglected: 



The Lagrange density can now be written aJ^ 

where r = t ■ cJo, oJq = STrM^/i^/fi, = uj/ujq. The dy- 
namics of SFi are then described by the following system 
of equations: 



(1 



N,, - TV, 



cos 2^ cos 2x + 26: I 7V^ - 7V3, + — I sin 2^ sin 2x 



(1 -e^) [{-iy - 7a,) COS 2^ + (7^ +73, ) COS 2^] - 4(/ia, + /i/ + e/i^) cos <l> cosx 



+ ^_e{hx -\- hf) -\- hd] sin<l>sinx — 4/i^ cosa;t(sin <l> cosx + £:cos<l>sinx] + 2^^^^ ^ + 2emzlz)- 



^-^(<l> + ex) = cos r^r(cos $ cosx — esin^sinx) — [hx -\- hf -\- ehd) sin $ cos x — [^{h^ ^ hf) ^ hd] cos $ sin x 



^{[(1 + e^) {Ny - N,) + hu] cos2x + (1 - e^){jy - 7,)} sin 2$ - e (^Ny - + 



COS 2$ sin2x, 



^2) 

(x + = ^t(£: cos <I> cos X — sin <l> sin x) ~ i^x h f -\- ehd) cos $ sin x — [^(^^c + ^/ ) + ^rf] sin $ cos x 

- ^ {[(1 + e') (A^^ - A^..) + /^n] cos 2$ - (1 - e^){-fy + 7^)} sin 2x - e (^Ny - + sin 2$ cos 2x, 



These equations naturally yield the two ground states 
in the absence of microwave field {hy =0), with <l> = 7r/2, 
X = i7r/2, corresponding to the antiferromagnetic vector 
values of the bilayer 1 = (±2, 0, 0). 

Solving Eqs. [S2] for small oscillations about the two 
ground states yields the frequencies of the collective spin 
resonances - the SFi eigen-modes: 

{n^f =Ny-Nx^jx^hu-hdsmxo±{{l-e^)-f^ 
+ [{hx + hf) sinxo - e{Ny - - 7:.)]'}'^'. 

Here "— " stands for acoustical (in-phase) and for op- 
tical (out-of-phase) resonance frequencies (Fig. [2|l). The 



criterion for generating spin resonance in both ground 
states, in terms of Hx, is given by 

(AttMs) 

I (l-£2)^2 

^ "^^ - 1 - hl-e^{Ny-Nx-jx)'' 

Here all terms are functions of (^1,^2) or (iYf^, i^2^). The 
optical resonance frequencies for the two ground states 
are plotted in fig. (2)3, c as surfaces in the SFi-asymmetry 
versus external field parameter space. The bottom 
plane of figs. (2)3, c are intensity plots of the difference 
in the optical resonance frequency of the two ground 
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FIG. 2: (Color online) a, Schematic of the acoustical (in- 
phase) and optical (out-of-phase) collective resonance modes; 
the horizontal dotted line represents the easy- axis; the colored 
shadows depict the individual macrospin oscillations. b,c, 
The splitting of the optical spin resonance frequency for the 
two AP states (color scale) as a function of: b, Hx and asym- 
metric biasing (fringing) field, with zero thickness imbalance 
; c, Hx and thickness imbalance, with zero fringing field. 



states (A/optical = Ar2+ • a;o/27r), again versus the key 
asymmetry- field parameters. 

In order to determine the mechanism of resonant SFi 
switching, we next analyze the behavior of the system for 



large microwave excitation amplitudes. For this we only 
consider a fringing field asymmetry (which typically has 
a stronger effect on the splitting of the optical resonance) 
and take the thicknesses of the two free layers to be equal 
{ti = t2 =^ e = 0). To simplify the derivations, the 
external uniform E A- field is set to hx = —hf, for which 
the resonance splitting is maximal. 

A hard- axis (y-axis) microwave field forces the two 
macrospins to oscillate out-of-phase, in the optical mode, 
which results in a modulation of x- At the same time, $ 
only slightly deviates from its ground state of <l>o = 7r/2. 
Therefore, it is a good approximation to take cp = ^—7t/2 
{\(p\ <C 1). As a result, Eqs. (S2) can be rewritten as 



2xa 



dx 
^ dr 



Q.^ sin X cos x 



(4) 



-^hd cos X + ^ cos VtT sin x = 0- 



) -^ = 0, 

2{Ny-N:, + K)-cos'^X 



(5) 



— /if^ • sin X + • cos VLt cos x, 



provided that Vt ^ VLq and ^ (which is 

valid in our case of elliptical particles of aspect ra- 
tio approximately 1). and VLq are the characteris- 
tic frequencies of small amplitude 'acoustical' and 'op- 
tical' oscillati ons in SAF, respectively, a nd are given 
by \/ {Ny - hu) - {^y - 7x) and Vt^ = 

Under these conditions. 



= ^{Ny -Nx^ hu) 

^{Ny-Nxihu)^{-{y^-{x)- 
the variation in can be regarded as 'slow' compared to 
the 'fast' oscillations in x, in which case averaging the x- 
and l]-terms in equation [S8] over one period is justified 
(represented by the bar). The Al^o-term in eq. S9 is the 
parameter characterizing the dissipation of energy in the 
system. Its value is associated with the Landau-Lifshitz- 
Gilbert damping parameter. 

Using Eqs. (S9) and (S8), the stability criteria for the 



ground states can be derived (see supplementary): 




where Ac is the critical oscillation amplitude derived 
in the supplementary, 7/ = (Q. — /^c^ with = 
- 3A2/4)/(l + 5A2/4). The above result takes 
into account that |7/| <C 1. 

Thus, in the presence of biasing field asymmetry, the 
stability regions for the two ground states split, and can 
fully separate such that microwave frequency-amplitude 
areas exist where only one states is stable, as illustrated 
in Fig. Is] If resonantly excited within these areas (shown 
in grey)7 the system switches into the given stable AP- 
state and remains in it, now off-resonance with the ex- 
citation. The separation of the minima in Fig. |3]is the 
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FIG. 3: (Color online) Stability diagram versus frequency- 
amplitude for the SFi ground states under resonant excita- 
tion. The curves separate the regions of stability and in- 
stability (below and above). The greyed out area shows the 
frequency-amplitude region where only one of the two states is 
stable, into which the system switches for given excitation. Q.o 
is the optical center frequency (that of an ideal SAP), hy the 
amplitude of the excitation, 7^ = (^2 — Qc)/^^c 7/ • 4Qo oc 
Q — Qc, ^ the excitation frequency. 



splitting in the optical eigen-mode of the system and is 
given by 7/ = hd/^l. 

The separation of the respective optical spin resonance 
frequencies is due to the difference in the effective biasing 
field acting on the individual layers. If the amplitude 
of the excitation is large enough and the frequency is 
chosen in resonance with a given ground state, switching 
into this state occurs, while the reverse switching is off- 
resonance and therefore disallowed energetically. 

Interestingly, although the oscillation of the two 
macrospins is of optical type (out-of-phase), the switch- 
ing occurs acoustically (in-phase rotation). The dipole- 
coupling is too strong for the moments to cross passed 
each other under the relatively weak microwave resonant 
fields used in this work. The switching trajectory is thus 
an in-phase rotation with superposed out-of-phase oscil- 
lations of smaller amplitude, with the energy effectively 
transferred from the pumped optical oscillations into a 
large- angle acoustical rotation. As a result, the expected 
fastest switching time is the half-period of the acousti- 
cal resonance frequency (corresponding to approximately 
500 ps for our experimental geometry). For still faster 
resonant switching one would aim in the SFi-device de- 
sign at maximizing the acoustical resonance frequency 
while maintaining a suitable size splitting in the opti- 
cal resonance frequency; for example, by increasing the 
aspect ratio of the nanomagnets. 



III. MACROSPIN SIMULATIONS 

The macrospin simulations in this section solve the LLG 
equations for the system in the time domain, taking into 
account all the external and internal forces acting on the 
magnetic moments of the two layers, without the small- 
signal assumptions used in obtaining the analytical re- 
sults. The simulation model is an extension of our pre- 
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FIG. 4: (Color online) a. Switching probability map for the 
sample parameters given in the figure, showing a clear separa- 
tion of the two AP states due to a biasing field and thickness 
asymmetry, b. Time trace of the two macrospins in the bilayer 
during a switching event, showing that the switching occurs 
via an in-phase rotation (acoustical mode), with superposed 
optical oscillations of smaller amplitude, c. Simulated easy- 
axis hysteresis for a sample with a thickness imbalance and 
biasing field asymmetry for s > (main panel), correspond- 
ing to our typical experimental configuration (left inset), and 
for £ < (top-right inset). In contrast, an ideal, symmet- 
ric spin-flop bilayer shows a strictly anti-symmetric response 
(bottom-right inset; simulated). 



vious work on the energetics of the SAF systerrP^. The 
magnetization vectors evolve in time driven by torques 
of the respective effective magnetic fields, with a suit- 
ably small time step of 5 ps. A thickness imbalance and 
fringing field are the two asymmetry parameters. The 
simulations include thermal agitation at T = 300 K. 

In the case of no thickness imbalance, with the only 
asymmetry due to the unequal biasing fields on the indi- 
vidual free layers (from the reference layer built-in into 
the nanopillar), the numerically simulated switching be- 
havior agrees well with that predicted analytically above. 
The optical resonance is clearly split, although the mag- 
nitude of the splitting is somewhat smaller. 
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A more pronounced resonance splitting is obtained by 
adding a thickness-imbalance to the magnetic asymme- 
try of the system. Figure [4^ shows a typical simulated 
switching map. The area where switching occurs only in 
one direction is where only one ground state is dynami- 
cally stable. The general form of the stability regions sim- 
ulated here for large-signal microwave excitation agrees 
well with the analytical results. 

Figure shows the real-time trace of the magneti- 
zation of the individual layers during a switching event. 
It is clear that the microwave-pumped oscillations are of 
optical nature but the switching itself is an "acoustical" 
in-phase rotation. The oscillations of the magnetization 
in the initial ground state under a resonant excitation 
of amplitude suitably large to produce switching are, in 
fact, smaller in the switched-to ground state since this 
state is now off-resonance with the microwave field. This 
effect of resonant state selection is a direct consequence 
of the magnetic asymmetry of SFi. 

The full macrospin model further yields important in- 
sights into the quasistatic magneto-resistance properties 
of the junctions. Using Eq. [Sljthe quasistatic switching 
(H^) and spin-flop fields {H2 ) can be obtained and used 
for extracting the asymmetry parameters (individual 
layers' thicknesses and effective biasing fields) from 
the measured data. The interrelation between the 
thicknesses and fringing fields is (in notations of Fig. 
gc)): 



--T{-me\{Ny-N,-^,)-hf 



± 



I I h+-h^ -ht + h2 -{h++h^ +4+^2) 

l''l'^4(iV,-iV,-7.)' 4 

Here = AnMshf is the field for which only one AP 
ground state becomes unstable and = AnMshf is 
the spin-flop field, at which both AP ground states be- 
come unstable. Note that the actual values of 
slightly smaller than the analytical values due to thermal 
fluctuations. 



IV. EXPERIMENTS 
A. Samples 

The samples used in the experiments were spin-flop 
type magnetic random access memory (TMRAM) cells 
with lateral dimensions of 450x375 and 420x350 nm, and 
elliptical in-plane shape. The stack consists of a magnet- 
ically soft SFi separated from a magnetically hard SAF 
reference layer by a thin Al-0 tunnel barrier. The spin- 
flop bilayer is composed of two dipole coupled permalloy 



ferromagnets (free layers) with thicknesses ti ^ t2 ^ ^ 
nm (excluding the magnetically 'dead' layers at the film 
surface^^ separated by a ~1 nm thin TaN spacer, for 
which there is no inter layer exchange coupling, J = 0. 
Here ti and ^2 are the effective thicknesses of the bottom 
and top free layer, respectively. The SFM reference layer 
is nearly ideally flux-closed and is therefore essentially 
insensitive to external fields of strengths used in the ex- 
periments discussed in this paper. The samples used in 
this work had a weak fringing field from the reference 
layer acting on the two free layers. The fringing field 
originated predominantly from the top fixed layer. The 
resistance of the stack was ^1 kQ, and the magnetoresis- 
tance -20%. 



B. Measurement setup 

In-plane quasi-static fields were applied using an exter- 
nal toroidal magnelP^. High-frequency fields were gen- 
erated by contacting the on-chip integrated 50 Q bit- 
and word-lines using surface probes with 0-40 GHz band- 
width. The bit- and word-lines were oriented at ± 45° 
with respect to the the easy-axis of the stack in the so- 
called "toggling" configuration. A current through the 
bit- or word-line results in an in-plane magnetic field 
perpendicular to the line. In the results presented be- 
low, the high frequency excitation was applied using the 
word line only, which was electrically decoupled from the 
stack. The applied high-frequency fields were AC pulses 
of desired frequency ranging from 1 MHz to 10 GHz. The 
AC fields were generated using an Agilent E8247C PSG 
CW Signal Generator with a bandwidth of 250 kHz - 20 
GHz. The wave-duration was controlled using an RF- 
switch driven by an Agilent 33250A AWG. The shortest 
possible pulse-length with this setup was 8 ns. The state 
of the sample was determined by measuring the resis- 
tance of the sample and thereby the relative orientation 
of the magnetization of the bottom free and top fixed 
layers. To separate the DC-signal from the RF-signal, a 
bias-tee was used. 



C. Experimental results 

The samples were SFi bilayer s of approximately 5 nm 
thick Permalloy, separated by 1 nm thick TaN, inte- 
grated into approximately 350x400 nm in-plane elliptical 
nanopillars, also containing an Al-0 read-out junction 
with a nominally flux-closed magnetic reference layer. 
The microwave excitation was applied using a closely 
spaced field-line making a 45° angle with the SFi easy 
axis. The ferromagnetic resonance frequencies were mea- 
sured using transport spectroscopy, where the junction 
resistance versus frequency showed a maximum or a min- 
imum for the low-resistance or high-resistance ground 
state, respectively. The resistance was measured while 
slowly sweeping the frequency of a continuous-wave ex- 
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FIG. 5: (Color online) a, Measured resonance frequencies of a 420x350 nm spin-flop bilayer. The intensity map (red and blue 
for the two ground states) is normalized to the maximum oscillation amplitude for the given EA-field. b, Switching probability 
for the two ground states subject to 300 ms long excitation pulses with -20 Oe static EA-field applied. Black color denotes 
the frequency- amplitude points, for which switching in both directions occurred. The inset shows the section of the analytical 
switching map (Fig. [sj covered in the measurement, c, Real time resistance traces of a sample excited by short microwave pulses 
of alternating frequencies in resonance with the split optical mode, showing controlled resonant state selection, d. Resonant 
state selection using resonant microwave pumping of fixed frequency and alternating external EA-field, tuning the two ground 
states in and out of the optical spin resonance. 



citation of fixed amplitude (typically 25 Oe). Figure [5|i 
shows the measured resonance frequencies for the two AP 
ground states of a SFi with Mq < 0, i^f ^ 7 Oe, ^ 3 
Oe and ti—t2 ~ 0.9 nm (values deduced from quasistatic 
EA- hysteresis as shown in figjl]^). The splitting of the op- 
tical spin resonance at ~ 4 GHz is in excellent agreement 
with the predicted behavior (Fig. |3|. 

Figure |5^ shows the measured oscillation amplitude as 
a function of the excitation frequency and applied bias 
field. The amplitude is normalized with respect to the 
maximum value for each bias field. This diagram clearly 
shows that both field-controlled and frequency-controlled 
switching is possible. 

Frequency-amplitude-switching maps were measured 
for both AP states (low- to high-resistance and high- 
to low-resistance) for different excitation pulse durations 
and static E A- fields. Figure [sJd shows the switching maps 
for both AP states of the soft bilayer under 300 ms pulses 
of resonant excitation, with = —20 Oe. Switching 
at frequencies lower than 3.9 GHz only occurs from the 
low- to high-resistance AP state {/l^h)^ while switch- 
ing at above 3.9 GHz only occurs from the high- to low- 
resistance state {/h^l)' Shortening the duration of the 



excitation pulse requires a higher amplitude for switch- 
ing the junction, at the same time leading to an increase 
in the error-rate likely due to a more non-linear and 
more non-uniform process for higher-amplitude pumping. 
Nevertheless, even for short excitations, the optical res- 
onance is clearly split and amplitude-frequency regions 
exist where the switching is strictly one-directional. The 
weaker higher-order oscillations in the probability maps 
of Fig. [sJd are due to secondary spin-wave excitations 
in our nanoparticles that are somewhat bigger than the 
true single-domain limit, which we also observe in our 
micromagnetic simulations as a small superposition on 
the well-defined macrospin behavior. 

Figure [5]:) shows a realtime trace of the junction resis- 
tance during the application of a sequence of microwave 
pulses of 100 ns duration, of frequency alternating be- 
tween the split optical spin-resonance peaks. A pair of 
pulses is sequentially applied at each frequency, demon- 
strating the strictly unidirectional character of resonant 
switching at a given frequency (for a given AP state of 
the SFi). 

Figure [sji shows a realtime resistance trace (blue) un- 
der continuous microwave excitation of fixed frequency 
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near the centre of the optical resonance, with the EA- 
field of amphtude sufficient for tuning the two AP ground 
states in and out of the resonance. Between the EA- 
pulses the excitation is turned off shortly to confirm the 
state the sample is in (microwave amplitude shown in red; 
EA- field in green). The EA- field pulses were generated 
with an external electromagnet, so their shortest dura- 
tion was 300 ms. The data show that the field-controlled 
resonance state selection is highly reliable. 

Both measurements were repeated many times with 
the resulting error-rate smaller than 10~^ (1000 mea- 
surements were performed without any error). The ef- 
fect is thus essentially deterministic for the sample and 
frequency-field parameters used, and still is expected to 
improve (e.g., in terms of speed) on scaling down the 
sample size. 

V. CONCLUSIONS 

We show how external microwave and static biasing fields 
can be used to tune the energetics and dynamics of a 
spin-flop system between symmetric and asymmetric be- 



havior. We demonstrate controlled state selection in a 
synthetic ferrimagnet using a frequency modulated reso- 
nant microwave field or a uniform external biasing field. 
This effect is explained analytically using a small-signal 
analysis of the SFi spin dynamics, predicting a tunable 
splitting of the optical spin resonance, which is in good 
agreement with our numerical arbitrary signal-strength 
simulations. The effect is fast, robust, and offers a new 
way of magnetic switching of spin-flop nanodevices used 
in such large scale applications as magnetic random ac- 
cess memory. 

Our analysis shows that the switching process is fast, 
^ 1/2 fa ^ 500 ps, and proceeds through an in-phase 
(acoustical) spin rotation with superposed out-of-phase 
(optical) spin oscillations. Experimentally we demon- 
strate essentially error-free controlled switching down to 
100 ns in pulse duration of the resonant microwave ex- 
citation. Smaller, more single- domain- like samples (sub- 
100 nm range), with optimized magnetic asymmetry and 
increased aspect ratio are expected to show significantly 
faster resonant switching due to reduced generation of 
unwanted spin-wave modes and higher acoustical fre- 
quencies. 
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Supplementary: Controlled resonant switching of a synthetic ferrimagnet 

I 



Macrospin theory 

Considering the free layers of the synthetic ferrimagnet 
(SFi) as single- domain particles, with the the z-axis out- 
of-plane, the EA and a uniaxial anisotropy field along 
the X-axis, the total magnetic energy {E) of the two free 
layers can be written affl: 

E ^ 

+ ljxmixm2x + ljymiym2y + Jjzmizm2z 
-2{hx + h'^)mix - 2hy cos {ujt)miy}, 

where Ms is the saturation magnetisation (assumed 
same in both layers); Vi = irabti/A the volume of 
layer i with a, b and ti the length, width and thick- 
ness of that layer, respectively; Nia the demagnetising 
factors of layer i in the a-dimension with Ni^ = 1 
and NjSx^yj = n^x,y}ti/b where n^^^yy are the reduced 
factors^ the cartesian unit magnetisation vector of 
the i-th layer; hu = Hu/^ttMs the normalised intrin- 
sic uniaxial anisotropy field along the easy-axis; hy = 
Hy/AirMs the normalised amplitude of the applied AC- 
field in the y-direction; u the frequency of the AC-field; 
hx = Hx/4:7tMs the normalised external easy-axis DC- 
field; = / AiiMs the normalised fringing field act- 
ing on layer i in the x-direction; 7^0, = r^Nicy the inter- 
layer exchange field from layer i acting on the other layer 
in the a-direction^. In our case ^io^ is solely of dipolar 



nature. Since in a thin magnetic film the in-plane de- 
magnetising factors are proportional to their thickness, 
they can be rewritten as 

Ni{xM = [1 - (-1)' ^] ^{^M^ = TV,, - Niy, 
where 

^ (ti + t2)n{a,^^} _ ti - t2 

The same can be done for the volume and interlayer ex- 
change field of the individual layers: 

Vi=[l- (-l)V] y/2; 7,^,,,} = [1 - [-Ife] 7{.,,}, 

For the description of the dynamics of the magnetisa- 
tion of the system, it is useful to convert the magnetisa- 
tion into polar coordinates: 

= {cosipi sin6>,, ^inipi sin6>,, cos6>,). 

Since in the used geometry the thicknesses are much 
smaller than the lateral dimensions and the sample is 
longer along the x-dimension than the y-dimension, the 
following holds: 

Nix < Niy < Ni, ^1; 6^= 7r/2 + ^i, l^il < 1, 

Hence, cos(^i) ^ — The terms r^, ^1, ^2 are minimal 
and therefore jzS,iS,j can be neglected: 



W-Wq 
7rM| 



+ TV, 



TV. 



cos 2<l> cos 2x^2e (Ny- N, 



sin 2<l> sin 2% 



(1-6^) [(7^ -7,)cos2c& + (7^ +7,)cos2<l>] - 4{h^ ^ hf ^ ehd) cos ^ cos x 



(SI) 



- 4:[e{hx -\- hf) -\- hd] sin $ sin x — ^hy cosa;t(sin $ cosx + ecos^sinx] + + ^2 + 2emzlz)^ 



where Wq = 7r[(l + e'^){Ny + TV,) + KjM'^; ^, x are the Hence, 
normal coordinates and are given by <l> = {(pi + ^2)/2, 
X = {^1 - ^2)/2; = ^1 + ^2, Iz = ii- 6, here 
and Iz are the normalised z-components of the magneti- 
sation vectors of the resulting pair of layers. Further 
the following relations were used: hf = (h'^' + /?'2^)/2, 
hd = (hT - h^)/2. 

The Lagrange density is then given hy^: 



\ dr dr 
dx d^ 



TiMs 



2^ 2MB 



(1 



cos^,)^(l 



where r = t ■ cjq, ^0 = ^'ttMs/j^b /^^ ft = uj/ujo. There- 
(— 1)V) — W. fore, the dynamics of the system are described by: 
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(p 

^-^($ + £X) = hy cosl]r(cos^cosx — esin^sinx) — {hx -\- hf -\- ehd) sin<l>cosx — [£{hx -\- hf) -\- hd] cos ^ sin x 

- ^{[(1 + s^) {Ny - Nx) + K] cos2x + (1 - e^){jy - 7,)} sin 2$ - e (^Ny - Nx ^ cos 2$ sin 2x, 

(P (S2) 
^-^(x + = hy COS r^r(£: COS <l> COS X — sin<l>sinx) — {hx -\- hf -\- ehd) cos<l>sinx — [£{hx -\- hf) -\- hd] sin<l>cosx 

- ^{[(1 + s^) {Ny - Nx) + hu] cos 2^ - (1 - e^){jy + 7,)} sin2x - e (^Ny - Nx ^ sin 2$ cos 2x, 



r 



From these equations it is also easy to determine that 
the ground states in the absence of the oscihating y- 
field {hy = 0) is given by <l> = 7r/2, x = i7r/2 which 
corresponds to the antiferromagnetic vector values 1 = 
(±2,0,0). 



criterion for Hx for degenerate resonance in both states 
is given by 



(1 - £2)^2 



/?2 



eHNy-N,--f,y 



A. Resonance frequency 

To switch the SFi resonantly between its two ground 
states, it is necessary to determine the resonance- 
frequencies of the ground states. If the amplitude of 
the oscillating field is small {\hy\ <C 1), then the mag- 
netic moments of the SFi deviate only minimal from their 
ground state: 

^ = 1 + ^1; x = xo + xi; xo = ±|; |^iUxi|<i 

(S3) 

Substituting the relations of (S3) into Eqs. ( |S2p gives: 



^22X1 



d^^i 
dr"^ 



- fli2Xi = shy COS Qt sin xo, 
^21^1 = cos fir sin Xo, 



(S4) 



ni = (1 + s'){Ny - Nx) + (1 - s'){jx + (-1)^,) 
+ [e{hx + hf) + hd] sinxo, 
= 2e{Ny -Nx^ hu/2) - {hx ^ hf ^ ehd) sinxo, 

where i 7^ j. The eigenfrequencies or resonance- 
frequencies of a system of equations as in Eqs. (S4) is 
easy to solve, and is given by 



{n^f = Ny-Nx^jx^hu-hd sinxo ± {( 1 - ) 7^ 
+ [{hx -^hf) sinxo - e{Ny - Nx - 7.)]'} 



2 
y 

1/2 



where "— " stands for acoustical (in-phase) and for 
optical (out-of-phase) resonance frequencies. Hence, the 



where all the terms are functions of (^1,^2) or (i^f^, i^2^). 
The difference in optical frequency of the two ground 
states (A/optical = Al^^ • ujo/27t) can be approximated 
by the linearisation of the frequencies around e = and 
hd = 0: 



hd±s{hx^hf)^^MZ^^ 



B. Resonant stability 

To determine the switching mechanism for the SFi us- 
ing resonant excitation we analyse the behaviour and sta- 
bility of the system for large excitation amplitudes. For 
this we only consider a fringing field asymmetry (which 
has the largest contribution to the splitting of resonance) 
and take the thicknesses of the two free layers to be equal 
(ti = t2 ^ £ = 0). To simplify the derivation, the ex- 
ternal DC field is set to hx = —hf (This is the external 
field at which the splitting is maximal). 

When an external AC-field of suitably high frequency 
is applied along the hard-axis (?/-axis) of the elliptical 
magnetic particles, the magnetic moments begin to oscil- 
late in the optical mode, which results in a modulation 
of X- At the same time, ^ only deviates slightly from 
its ground state $0 = 7r/2. Therefore, it is a good ap- 
proximation to take (p = ^ — it/2 {\lp\ <C 1). Taking all 
the above assumptions in consideration, Eqs. (S2) can 
be rewritten as 
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[nl - 2{Ny -N, + hu) cos\ 

— hd sinx + hy cos Qt cosx] f = 0, 



(S5) 



sinx cosx + cosx + ^2/ cos ftr sinx = 0, (S6) 



with and the characteristic frequency of smah 
amphtude 'acoustical' and 'optical' oscillations, re- 
spectively, in a symmetric system and are given 
by = \/ {Ny - hu) - {jy - Jx) and = 

y^{Ny -N^i hu) + {jy + -fx)- Provided that l^a < ^o, 
(which is valid for elliptical particles with aspect ratio 
of approximately 1) the variation in cp can be regarded 
as 'slow' compared to the 'fast' oscillations of x- If the 
frequency of the applied external field is near the optical 
resonance frequency l^o, then averaging the X" and Q- 
terms in Eq. (S5) over one period is justified and yields 



{nf)'-^ = o, 



ni-2{Ny-N,^hu)-cos^X 
— hd - sinx -\- hy ■ cos cos x, 



(S7) 



(S8) 



where the bar stands for averaging over one period of 
the 'fast' oscillation. Using Eq. (S7) the stability of dif- 



ferent magnetic configurations resulting from the action 
of the AC-field excitation can be determined. Whenever 
fl^ff becomes complex, the state of the particle becomes 
unstable, which can result in switching. 

To evaluate Eq. (S6) for external excitations at the 



optical frequency we add a phenomenological term de- 
scribing the dissipative processes in the system: 



o sm X cos X 
cos X + cosier sinx = 0. (S9) 



Here XQq is the parameter characterising the dissipation 
of energy in the system. Its value can be associated with 
a dissipative Gilbert parameter. 

Starting in the ground states = 7r/2 and x = i7r/2) 
and turning on the AC-field results in periodic oscilla- 
tions of x, which will be described using a new variable 
/, connected with x by the following relation: 



X = Xo + 2 arctan(/), xo = ±7r/2 



(SIO) 



Here and further on the upper sign corresponds to the 
state in which the magnetic moments fiuctuate in the 
vicinity of / = (2, 0, 0), and the lower sign to fiuctuations 
in the vicinity of / = (—2, 0, 0). 



Substituting Eq. (|S1Q[) into Eq. (IS9| results in a fourth 



order differential equation: 



{l + f)(-^+2Xn,^+nl-hasinxo]f 



/+^sinxocosl^r(l-/4) = 0. 



(Sll) 



To reduce the degree of nonlinear terms in Eq. (Sll), 
the accuracy of the theory is limited to 

f « 1. 



This requirement is not too restrictive and allows to con- 
sider nonlinear oscillations with high amplitude. Follow- 
ing the standard procedure for finding the resonance in 
a nonlinear system, such as Eq. (Sll), the solution is 
represented in the form 



/ = Acos{nT^S)^C, 



(S12) 



with A the amplitude of the main harmonic, 6 the phase- 
shift with respect to the applied field, and ( an additional 
term taking into account anharmonic corrections. 

When the external excitation frequency Q is close to 
the resonance frequency, then the main harmonic oscil- 
lations dominate and \(\ <C A, in which case the an- 
harmonic corrections can be neglected. Substituting Eq. 
(S12) into Eq. (Sll ) and only considering the terms pro- 



portional to the main harmonics, (sin {^r) and cos (^^r)), 
the resulting amplitude and phase-shift of the oscillations 
become: 



3_ 

4^2 



Pl-hd sinxo) (1 



3A^ 



tan (5 = 



9/ 

(2A!^!^o) (1 + ^ 



(S13) 



i^2(l + 5^2/4)_(f22. 



-/idsinxo)(l-3A2/4)" 

(S14) 

From these equations it follows that at the excitation fre- 
quency ofn = y/{ni-hd sin Xo ) ( 1 - 3 A74) / ( 1 + 5 /4) 
the system is in resonance when the dissipation is weak. 
In this case the amplitude reaches a maximum Amax de- 
fined by Amax{l + AI^^JA) = hy/AXQQo- 

Substituting Eq. (|S12D into (|S1Q|) into ((SSl) and tak- 



ing the appropriate averages gives the effective 'acoustic' 
resonance frequency and with that the stability of the 
ground states: 
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{Ny - iV, + hu)4A^ 

(1+^2)3/2 



hd sinxo 



2 cos (5 sin ^1 



1 



(S15) 



r 



When the oscihation amphtude A substantiahy increases, 
the right hand side of Eq. (S15 ) becomes negative, which 



makes the respective state unstable. The difference in the 
resonance frequency for the two ground states / = (2, 0, 0) 
and / = (—2,0,0) is proportional to the asymmetry in 
the fringing field (hd). Therefore, it is possible to deter- 
mine the region in the frequency- amplitude phase space, 
in which one state is unstable while the other is stable. 
It is this asymmetry-induced selectivity in, for example, 
frequency for a given amplitude that allows to uniquely 
select a given AP state of the SFi. 

To determine the frequency-amplitude range where 
only one of the states is stable, the critical amplitude 
at which ft^^ becomes imaginary needs to be determined 
for the the two AP states. To do this, Eq. (S15) needs to 
be solved for the case where = 0. Assuming that the 
fringing field difference and the AC-fi eld amplitude are 
small, then the last two terms in (S15) can be neglected. 



and the critical amplitude Ac becomes: 



4Al {l^Al 



-3/2 ^ {Ny - N^^ K) - (7?y - Ix) 

{Ny - AT, + hu) 



Here the out-of-plane terms are neglected since we are 
only interested in the in plane AP ground states in which 
case Eq. (S16) is the resulting equation. The two ground 

7r/2 and x 



states are given by <l> = <l>o = tt/ 'Z an d X = Xo 
In this case the solutions of Eq. (S16) are: 



±7r/2. 



hx = e{Ny - Nx- jx) sin xo - hf 



± 



{Ny -N^^j^^K-hd sin xo)' - (1 - e^)-/^ 



(Sir) 

Now since the dipole coupling is strong and the aspect 
ratio relatively close to 1, Ny — N^ — jx < Therefore 
the sign of e determines which of the states have positive 
or negative spin-fiop or switch fields. Hence, in order 
to accurately give the relation for t he different transi- 
tion fields {H^ ^ H^)^ equation (S17) has to be slightly 
altered: 



AttMs 
± 



T{-ir\e\{Ny-N,-^,)-hf 



(S18) 



which reduces to a cubic equation and is straightforward 
to solve. For our samples the equation has three real 
solutions, of which one, A^ = 0.118, satisfies the criterion 
of applicability, <C 1. 

Once Ac is determined, Eqs. ( |S13D and ( |S14D can and fringing fields, Eq. (|S18|) has to be rewritten in the 



Now to be able to determine the individual thicknesses 



be used to determine the frequency-amplitude values for 
which each state changes between stable and unstable. 
Hence, a ground state is only stable if 



fundamental parameters, such that none are dependent 
on thicknesses or fringing fields. This results in the fol- 
lowing set of equations: 



hy_ 



where 7j 




7/ + 



hd sinxo 



+A2 1 



AttMs 



= T(-1) 



{ny-{l+rx)nx)- 



.3A2/4)/(l + 5^2/4). 
into account that |7/| <C 1. 



with the critical frequency Qc = 
The above result takes 



± 



(tl+t2) 



2b 

{ny-{l-rx)nx) 
2b 



SttMs 



±{-iY 



C. Switch and spin-flop fields 



ll — l2 -til ~ ^2 

\h-t2\ SttMs 



'Mm 



.1/2 



The switching and spin-flop flelds that are determined 
by measuring the EA hysteresis of the sample are used to 
calculate the approximate thicknesses and fringing-fields. 
The derivation of this formula is quite straightforward. 
The switching and spin-fiop occur when the initial ground 
state is no longer an energy-minimum. Starting from the 
energy-density given by Eq. ( SI ), the switching and spin- 



fiop fields can be determined by solving hx for 



d^dx J 



(S16) 



t2\ 



{uy - (1 + ra;)nx) SttMs ' 



■H^)/2. 
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